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ABSTRACT 

The thin red giant branch (RGB) of the Carina dwarf spheroidal galaxy appears at first sight quite puzzling 
and seemingly in contrast with the presence of several distinct bursts of star formation. In this Letter, we provide 
a measurement of the color spread of red giant stars in Carina based on new BVI wide-field observations, and 
model the width of the RGB by means of synthetic color-magnitude diagrams. The measured color spread, 
(jy-i = 0.021 ±0.005, is quite naturally accounted for by the star-formation history of the galaxy. The thin RGB 
appears to be essentially related to the limited age range of its dominant stellar populations, with no need for 
a metallicity dispersion at a given age. This result is relatively robust with respect to changes in the assumed 
age-metallicity relation, as long as the mean metallicity over the galaxy lifetime matches the observed value 
([Fe/H]= -L91 ±0.12 after correction for the age effects). This analysis of photometric data also sets some 
constraints on the chemical evolution of Carina by indicating that the chemical abundance of the interstellar 
medium in Carina remained low throughout each episode of star formation even though these episodes occurred 
over many Gyr 

Subject headings: galaxies: individual (Carina) — galaxies: dwarf — galaxies: evolution — galaxies: stellar 
content — Local Group 



in particular the role of galaxy-wide gas outflows and/or infall, 
remains a controversial issue and many basic questions are still 
open. Recent efforts to model the chemical evolution of dwarf 
spheroidal galaxies have been presented by Ikuta & Arimoto 
(2002), Carigi, Hernandez, & Gilmore (2002), and Tolstoy et 
al. (2003). All groups use known SFHs to infer the chemical- 
enrichment histories of dSph galaxies, but the scenarios differ 
in the role played by galactic winds and gas infall. 

This Letter shows that realistic modeling of the color distri- 
bution of red giant stars by population synthesis methods can 
help constraining the chemical evolution of dwarf spheroidals 
galaxies. The intrinsic color width of the upper RGB in Ca- 
rina is quantified here by exploiting the high statistics of a new 
BVI database from wide-field observations. The star-formation 
history of Carina is then derived for a set of different metal- 
enrichment histories (the age-metallicity relation of this galaxy 
is still uncertain) by comparing observed and synthetic color- 
magnitude diagrams. The narrow RGB of Carina is shown to be 
consistent with the star-formation history derived from the ob- 
served CMDs, under reasonable assumptions on the chemical- 
enrichment history. 

2. THE DATA 

Observations of Carina were carried out using the Wide Field 
Imager at the ESO/MPG 2.2m telescope. The camera consists 
of a mosaic of eight 2k x 4k EEV CCDs, in a square array of 2 
X 4 detectors. In this configuration, the total number of pixels 
is 8k X 8k, with a scale of 0.238 arcsec pixel"' and a field of 
view 32' X 32'. Three different fields were observed: a central 
field and two outer fields with ~ 20' offsets. Sky conditions 
were photometric, and standard star fields taken from the list of 

' Based on data collected at E.S.O. La Silla, Chile, Proposals No. 63.N-0017 and 64.N-0512 

^ INAF, Osservatorio Astronomico di Padova, vicolo dell'Osservatorio 5, 1-35122 Padova, Italy; (held,i'izzi)@pd.astro.it 
^ Dipartimento di Astronomia, Universita di Padova, vicolo dell'Osservatorio 2, 1-35122 Padova, Italy 

Istituto di Astrofisica Spaziale e Fisica Cosmica, CNR, Area di Ricerca Tor Vergata, Via del Fosso del Cavaliere 100, 00133 Roma, Italy 
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L INTRODUCTION 

Since its discovery in 1977 (Cannon, Hawarder, & Tritton 
1977), the Carina dwarf spheroidal (dSph) galaxy has played 
an important role in our understanding of the nature of dwarf 
galaxies. Carina was in fact the first dSph galaxy in which a 
prominent intermediate age population was recognized (Mould 
& Aaronson 1983). The presence of a non-negligible compo- 
nent of old stars was subsequently revealed by the detection 
of a large number of RR Lyrae variable stars (Saha, Monet, & 
Seitzer 1986; DaU'Ora et al. 2003). Mighell (1990) estimated 
the old population to constitute (17 ±4)% of the galaxy popu- 
lation, a number confirmed by Smecker-Hane et al. (1994) and 
Mighell (1997). A detailed description of the star-formation 
history (SFH) of Carina was given by Smecker-Hane et al. 
(1994, 1996). A deep BR color-magnitude diagram (CMD) re- 
vealed the presence of three distinct turn-offs, pointing at three 
separate episodes of star formation approximately 2,3-6 and 
11-13 Gyr ago (see also Monelli et al. 2003). A quantitative 
assessment of this scenario was then provided by CMD simu- 
lation techniques applied to CTIO 4m and HST data (Hurley- 
KeUer, Mateo, & Nemec 1998; Hernandez, Gilmore, & Valls- 
Gabaud 2000; Dolphin 2002). 

Quite surprisingly given this complex star-formation history. 
Carina shows a very thin red giant branch with all the differ- 
ent turn-offs connecting at the RGB base (Smecker-Hane et 
al. 1994, 1996). This seems to imply that either no significant 
metal enrichment took place in Carina, or that age and chemical 
enrichment of the stellar populations played in opposite direc- 
tions to yield a narrow RGB. While these alternative views have 
often been discussed, they have never been quantitatively ad- 
dressed. The chemical evolution of dwarf spheroidal galaxies. 
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Landolt (1992) were observed for calibration purpose. A full 
description of the data and reduction methods will be presented 
in a forthcoming paper (E. V. Held et al. 2003, in prep.) while 
just a brief outline is given here. Pre-reduction of CCD images 
was performed within the IRAF environment. Each image was 
bias- subtracted and flat-fielded using twiUght sky flats. After 
these steps, aU images were astrometrically calibrated using the 
IRAF MSCRED package (Valdes 1998) and the script package 
WFPRED developed at the Padova Observatory. 

Photometry was obtained using DAOPHOT Il/ALLSTAR (Stet- 
son 1987) on the co-added images, and the results were cali- 
brated using calibration relations computed from standard stars 
observed in all CCDs. Similarly, aperture corrections were 
computed using a growth curve analysis performed on bright 
isolated stars in each CCD. Artificial star tests were done to 
evaluate the uncertainty of photometry and the degree of com- 
pleteness of the data. The photometric errors in the V band 
range from ^ 0.01 mag near the tip of the RGB, to ~ 0.3 mag 
near the detection limit. 

3. MODELING THE RED GIANT BRANCH OF CARINA 

3.1. RGB width measurement 

These wide field data allowed us to obtain a sound quantita- 
tive estimate of the RGB width in Carina. In order to measure 
the color spread of the RGB, a fiducial ridge line was fitted 
to the mean colors of red giants in 0.5 mag bins. The distribu- 
tion of color residuals about the fiducial line was then measured 
down to 1.5 mag below the tip. After quadratic subtraction of 
the contribution of photometric errors, the estimated intrinsic 
color spread of the upper RGB is: 

(7(v-/)= 0.021 ±0.005 

where the statistical error was estimated by bootstrap resam- 
pling of the RGB sample in 500 trials. 

Using a direct comparison of the red giant branch of Carina 
with the ridge lines of Galactic globular clusters (Da Costa & 
Armandroff 1990), we estimated a mean metalhcity [Fe/H] 
-2.08 ± 0.12 on the Zinn & West (1984) scale (E. V. Held et al. 
2003, in prep.). This value agrees well with previous photomet- 
ric determinations (Mould & Aaronson 1983; Mighell 1990; 
Smecker-Hane et al. 1994; Mateo, Hurley-KeUer, & Nemec 
1998), as well as the low dispersion spectroscopic measure- 
ments of the Call lines of Da Costa (1994) and Smecker-Hane 
et al. (1999). A somewhat higher metallicity, [Fe/H] = -1.6 
with a spread of 0.5 dex, has recently been derived by Tolstoy 
et al. (2003) from high-resolution spectroscopy of 5 red giant 
stars. 

The difference between this new spectroscopic estimate and 
previous photometric [Fe/H] value can be only partly accounted 
for by age effects on the RGB color. The stars in Carina are, in 
fact, ^ 7 Gyr younger than the stellar populations of Galac- 
tic globular clusters (see below), and we know that younger 
isochrones are shifted blueward (see, e.g.. Da Costa 1998). Us- 
ing Girardi's et al. (2000) models in the range Z = 0.0001 to 
0.0004, we find the effects of a younger age on the (V-/) color 
to be ~ 0.035 mag at the RGB level used for metallicity esti- 
mation. This bluer color mimics a lower metallicity, by ~ 0.17 
dex (see also Smecker-Hane et al. 1994). By taking this ef- 
fect into account, the age-corrected metaUicity of Carina is then 
[Fe/H]=-1.91±0.12. 




(v-i) (v-i) 

Fig. 1. — Comparison of the observed and synthetic color-magnitude dia- 
grams for Carina. The left panel shows the statistically decontaminated CMD 
of the central region of Carina, along with the boxes used by the minimiza- 
tion technique. The right panel shows the simulated color-magnitude diagram, 
assuming the PT law described in the text. 

The color distribution of RGB stars formally corresponds 
(for an old population) to a metallicity dispersion a-pe/H] = 
0.10 ± 0.04, which is significantly smaller than the quoted dis- 
persion of spectroscopic measurements. 



3.2. Population synthesis 

Figure 1 (left panel) shows the CMD of the the central part 
of Carina, for which simulations were performed. An elUpti- 
cal region with a 12 arcmin semi-major axis, 65° position an- 
gle, and 0.33 ellipticity (Irwin & Hatzidimitriou 1995) was se- 
lected to this purpose. The intrinsically poor statistics of the 
red giant stars in Carina prevents reliable measurements of the 
RGB width in the outer regions. The observed CMD was simu- 
lated using a parametric technique similar to that employed by 
Aparicio, Gallart, & Bertelli (1997), and a description of our 
methods is presented by Rizzi et al. (2002). In brief, our tech- 
nique assumes a number of known input parameters, such as 
the initial mass function (IMF), the distance to the galaxy, and 
the reddening. A Salpeter (1955) IMF is adopted, since recent 
work suggests a general validity of this IMF for masses larger 
than 0.5Mq (Scalo 1998; Kroupa 2000). A distance modulus 
(ot-M)o = 20.00 ±0.06 and a reddening £(B-y) = 0.045 ±0.01 
were adopted (E. V. Held et al. 2003, in prep.). Using these in- 
put parameters, a set of synthetic stellar populations were gen- 
erated using models from the library of Girardi et al. (2000) 
and the simulation code ZVAR that linearly interpolates the 
evolutionary tracks both in age and metallicity, including the 
effects of binary stars. The observational effects were applied 
to the theoretical random populations using the photometric er- 
rors and completeness obtained from the results of artificial star 
experiments. 
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Fig. 2. — The different chemical-enrichment histories tested in this work. 
The continuous line represents the adopted PT law, while the long-dashed and 
short-dashed lines show the linear and constant laws, respectively. 

A basic ingredient of the simulation is the chemical- 
enrichment history (CEH), for which we chose to test a num- 
ber of reasonable assumptions (Fig. 2). The first tested law is 
characterized by a prompt initial enrichment at ages of ^ 12 
Gyr, which is responsible for producing much of the observed 
metals, followed by a slow metallicity increase leading to a fi- 
nal value of [Fe/H] ~ -1.9 (hereafter "PT law", after the age- 
metallicity relation modeled by Pagel & Tautvaisiene 1998 for 
the Large Magellanic Cloud). This trend is similar to the 
CEH recently derived from spectroscopic observations of dwarf 
spheroidal galaxies (Gallart et al. 2002; Tolstoy et al. 2003). 
The final metallicity was constrained to yield a mean metal- 
licity over the galaxy lifetime equal to the value derived from 
the photometry of Carina. This was done iteratively since the 
galaxy's star-formation history is affected by the CEH normal- 
ization, and vice versa. 

Other tested laws, illustrated in Fig. 2, include a constant 
metallicity law at [Fe/H] = -1.91 and a hnear law starting 
at [Fe/H] = -2.3 and reaching [Fe/H] = -1.76 at the present 
epoch. These laws sample a wide range of possibilities in 
the rate of metal enrichment, still yielding the correct mean 
[Fe/H] over the SFH. Other chemical evolution laws with a 
more metal-rich normalization would be inconsistent with the 
photometric data. 

A total of 7 populations were generated using these input 
laws, with age bins centered at 1, 2.25, 4, 6, 8.5, 1 1 and 13 Gyr, 
each bin containing 30.000 stars. Tests were made with both 
a 30% fraction of binary stars and no binaries at all. The SFH 
that best fits the observed color-magnitude diagram was found 
by comparing the number of real and simulated stars inside a 
number of regions in the CMD ("boxes"). These were chosen 
to map (i) the luminosity function along the main sequence, (ii) 
the relative contribution of blue and red horizontal branch stars, 
(iii) the number of "red clump" stars, and (iv) the total number 
of stars along the RGB (see Fig. 1 , left panel). The width of the 
RGB itself is not an input constraint for the best-fit CMD. 



Fig. 3. — Relative star-formation history for the central region of Carina. 
The rates are normalized to a lifetime average rate of 25.6 Mq Myr"' . The er- 
ror bars represent 3(T confidence intervals from the chi-square statistic, while 
dashed lines show the upper and lower limits to the SFH derived from 100 
repeat CMD simulations. The bulk of the star formation appears to have hap- 
pened in a long episode between ~ 3 and ~ 7 Gyr ago. 

The best fitting synthetic diagram, obtained by minimizing 
the differences in star counts between the simulated and ob- 
served CMD with a downhill minimization algorithm, is 
shown in the right panel of Fig. 1 for the PT law. 

The derived star-formation histories are shown in Fig. 3, for 
the different assumptions on the metal-enrichment history. Our 
analysis of wide-field data agrees with previous results in that 
^ 75% of the stars formed between 3 and 7 Gyr ago, with a 
small yet significant old stellar component. Note that, although 
the details of the SFH are different, the main episodes of star 
formation are relatively independent of the adopted CEH. In 
fact, our suitable choice of the CMD regions over which stars 
are counted makes the derived SFH more sensitive to the life- 
times of evolutionary phases than to the metal enrichment law 
(and uncertainties on the isochrone colors). 

3.3. A thin RGB in Carina 

We can now return to the problem of the origin of the ex- 
tremely thin RGB in Carina, having in hand the necessary 
tools to quantitatively evaluate the effects of the chemical- 
enrichment and star-formation history. The color distribution 
of simulated RGB stars about a fiducial line was measured for 
each choice of the CEH, using the same methods as for the real 
data. We did not make use of the absolute colors of the syn- 
thetic RGBs, because of known uncertainties on the isochrone 
colors. The mean widths measured for the simulated diagrams 
are (T(v_/) = 0.023 ± 0.003 for the linear and constant laws, and 
a(v-i) = 0.025 ± 0.003 for the PT law. The mean values and 
statistical errors result from 300 trials. These results are es- 
sentially independent of the adopted fraction of binary stars. 
The simulations thus indicate that the Carina's narrow RGB 
naturally arises from the SFH derived from the data. The re- 
sults for all the tested enrichment laws agree with the observed 
color spread within one standard deviation. Indeed, these val- 
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ues probably mark the lower limit to the measurable RGB width 
set by stochastic effects and photometric errors in the CMD 
simulations. As expected, no conclusions can be reached on the 
details of chemical enrichment from photometry alone. A large 
sample of spectroscopic measurements will help, provided that 
it is of the highest precision and accuracy. 

4. DISCUSSION AND CONCLUSIONS 

In this Letter, we have reported new measurements of the 
mean metalUcity ([Fe/H]=-1.91 ±0.12, after correction for the 
effects of age) and the color spread (uv-/ = 0.021 ± 0.005) of 
red giant stars in Carina, based on new wide-field observations. 
Using synthetic color-magnitude diagrams to model the color 
distribution of RGB stars, we have quantitatively shown that 
the thinness of the red giant branch is a direct consequence of 
the derived star-formation history, for most chemical evolution 
laws that are plausible and consistent with the observed mean 
metallicity of the galaxy. 

The most fundamental reason for the RGB being so thin is 
the fact that the predominant intermediate-age stellar popula- 
tion that makes up the RGB, was formed in a limited time in- 
terval (5 ± 2 Gyr ago), with only a minor contribution by an old 
population. In theory, both the age range and the metallicity 
evolution of the red giant stars contribute in some way to the 
color range of the RGB, but the effects are subtle and hardly 
measurable. In fact, (i) the age interval is too small for the 
isochrone separation to be measurable, and (ii) the sensitivity 
of the RGB color to changes in [Fe/H] is low for metal-poor 
systems, so the effect of modest metal enrichment during the 
main star-formation episode is negligible too. 

Our photometric analysis bears some interesting implications 
on the chemical evolution history, by providing evidence that 
Carina could not experience a significant metal enrichment (i.e. 



above [Fe/H]~ -1-9), certainly not before 3 Gyr ago. In fact, 
the lifetime-averaged metallicity inferred from photometric ob- 
servations definitely rules out a metal-rich final composition. 
Little information is available on the most recent epoch, be- 
cause stars younger than ^ 3 Gyr are only a trace population, 
but for the same reason no significant metal enrichment is ex- 
pected in the last few Gyr. The most likely scenario remains that 
indicated by chemical evolution models, in which substantial 
metal enrichment is associated with the earliest star-formation 
epoch (> 12 Gyr). 

In summary, both the episodic star-formation activity (with 
long-lasting episodes and quiescent phases) and the lack of a 
significant metal enrichment suggest that metal-enriched gas 
outflows play a role in the evolution of Carina and other low- 
mass galaxies (Vader 1986; Mac Low & Ferrara 1999; see also 
Gallagher & Wyse 1994; Kunfli & Osthn 2000; Carigi et al. 
2002). Closed-box models (e.g., Ikuta & Arimoto 2002) ap- 
pear unable to explain the very low metal abundances of dwarf 
galaxies like Carina. In a similar way, the oxygen abundances 
measured in the interstellar gas of dwarf irregular galaxies ap- 
pear hardly consistent with a closed-box model (see, e.g., Sa- 
viane et al. 2002, and references therein). Observational ev- 
idence of galactic winds provides support to this scenario for 
active starbursts, a recent example being the Chandra obser- 
vation of a metal-enriched galactic wind in the dwarf galaxy 
NGC 1569 (Martin, Kobuhiicky, & Heckman 2002). 

We thank the referee for suggesting some significant im- 
provements to the paper, Drs. L. Carigi and X. Hemadez for 
helpful remarks, and M. Tosi for useful conversations on dwarf 
galaxy evolution. Support for this work was provided by the 
Itahan Ministry for University and Research through grants 
COFIN2001028897 and COFIN2002028935. 
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